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In this letter we show how to perform a systematic perturbative approach for the mode-coupling 
theory. The results coincide with those obtained via the replica approach. The upper critical 
dimension turns out to be always 8 and the correlations have a double pole in momentum space 
in perturbations theory. Non-perturbative effects are found to be very important. We suggest a 
possible framework to compute these effects. 
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In the mean field theory of glasses there are two transi- 
tions, the dynamical-replica transition, that corresponds 
to the mode-coupling transition in the dynamics and hap- 
pens at a temperature T c , and a thermodynamical transi- 
tion (the Kauzmann transition) that happens at a lower 
temperature Tk [H 0] ■ The dynamical- replica transition 
can be identified by looking at equilibrium properties of 
the system, i.e. its landscape; it corresponds to the for- 
mation of local minima in the free energy and it is usually 
studied using replicas [3-5]. The mode-coupling transi- 
tion is defined by the dynamical properties of the system. 
The two transitions are related as far as in the mean field 
approximation the time needed to escape from a local 
minimum of the free energy is infinite. 

This picture is exact in many solvable models. However 
it should be modified in the real world, where the mean 
field approximation is no more exact. In this letter we do 
not address the fate of the thermodynamical transition 
and we concentrate our attention on the dynamical mode- 
coupling transition. 

It is quite evident that in short range systems there 
is no dynamical transition, exactly for the same rea- 
sons for which there are no infinite lifetime metastable 
states. However if we neglect the so-called "activated" 
process, the dynamical transition is present. Moreover 
there is a very large amount of experimental and numer- 
ical data that are well fitted by the predictions of the 
mode-coupling theory, so that it is certain interesting to 
try to understand which is the critical behaviour associ- 
ated to the mode-coupling transition. 

In this letter we present a computation of the upper 
critical dimensions and of the critical properties of the 
dynamical mode-coupling transition: in the dynamics we 
consider only the mutual dependence of quantities that 
do not depend explicitly on the time, i.e. time has been 
eliminated parametrically as it happens in the general- 
ized fluctuation dissipation relations [(| @] • We will firstly 
present the results in the framework of the equilibrium 
replica approach and we will later show how the same 
results hold for the mode-coupling transition. 

The critical behaviour at the dynamical transition 
stems from the presence of dynamical heterogeneities [g- 



13]. In the dynamics these heterogeneities are related 
to the presence of correlated movements of cooperatively 
rearranging regions [l4| that have been observed both 
above and below the critical region around T c . We are 
interested in getting precise predictions on the properties 
of dynamical heterogeneities. 

Let us start with the basic definitions. Given two con- 
figurations of the coordinates (that we label with a and 
r), we indicate with q a ^ T (x) the similarity (overlap) of 
the two configurations in the region of space around the 
point x (many different definitions are possible). Usually 
q is equal to one for identical configurations and it takes 
a small value for uncorrelated configurations 
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For example we can take q (7>T (x) to be one if a region 
around x of size a _3.li] has the same particle content in 
the configurations a and r; otherwise q a ^ T (x) = 0, if the 
particle content is different. 

Let us consider the case where a is an equilibrium con- 
figuration of the system and r(t) is a configuration ob- 
tained using some dynamics at time t starting from the 
a configuration (i.e. r(0) = a). If the dynamics is non- 
deterministic, the configuration r(t) will depend also on 
some extra random variables rj. For simplicity of notation 
we will not indicate the dependence of r(t) on 77, unless 
we need it in an explicit way. We can define 



C(t) = q a (x,t) where q a {x,t) = {q a ,r(t){x)) ■ (1) 

Here the overline denotes the average over the Boltzmann 
distribution of the initial configuration (a) and the angu- 
lar brackets the average over r\. C(t) is the usual equal 
point (smeared over a region of size a) density-density 
correlation. Approaching the dynamical transition, C(t) 
will decay slower and slower, and will also develop a 
plateaux (as function of ln(i)) at the value Cp. This 
plateaux becomes infinitely long at the mode-coupling 
temperature; below the mode-coupling temperature, ne- 
glecting activated processes, the correlation does not de- 
cay any more, i.e. lim t _ i . 00 C(t) = C^ > Cp > 0. 

For the study of dynamical heterogeneities it is usual 
to consider a dynamical susceptibility Xi(~k) defined as 



V X i{t) = Q a (t)* - Q a {t) = (Q a (t) - C(t)f 



(2) 
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where V is the volume of the system and Q a (t) is the 
space integral of q<r{x, t). The quantity Xi{t) 1S a measure 
of the differences that are observed during the evolution 
from region to region In a similar fashion we can 

define the time dependent correlation function 



Giix -y,t) = q a (x, t)q a {y, t) - C(ty 



(3) 



The explicit time dependence is quite a complex prob- 
lem that we do not address in this letter; different dy- 
namical exponents are involved and they are known to 
be not universal (at least in the mean field approxima- 
tion) . We will consider here only relations where the time 
is not explicitly presents, as the dependence of Xi on Cj 
that can be obtained by eliminating the time paramet- 
rically, e.g. by plotting x<i(*) versus C(t) exactly in the 
same way as for the fluctuation dissipation relations. A 
particular example of a time independent quantity is x\i 
that is defined as the maximum of Xi(^) (that happens 
at time t* ) (the corresponding correlation will be denoted 
by G%(x)). Analogously we define a C-dependent corre- 
lation function G^{x — y\C) as follows 



Gi(x - y\C{t)) = G 4 {x - y,t) 



(4) 



We are interested in the universal properties of the pre- 
vious quantities when we approach the dynamical tran- 
sition. We denote by e the distance in temperature from 
the dynamical transition and we indicate from here on 
the dependence on e. We are interested in the behaviour 
in the double scaling limit e — >• and C —> Cp. Within 
mean field theory, the susceptibility X4(Cp, e) is divergent 
when e — > from above [32j . Given the obvious relations 
Xt( e ) = I Gt( x > e ) dx the divergence of xti e ) implies the 
existence of a divergent correlation length £(e). In the 
same way as in standard phase transitions, we expect the 
following scaling laws 



Xl(t) « 



Gl(k,e) = xl(e)Hk^)) , 



(5) 

where G|(fc, e) is the Fourier transform of G%{x, e). 

What is the replica counterpart of this behaviour? In 
the replica approach [H| we consider two replicas a and r. 
We denote by H(a) the original Hamiltonian, while the 
Hamiltonian of the r system is 



H(t) - hQ a 



(6) 



The thermal averages are taken first with respect to the 
t variables and later with respect to a variables. The 
dynamical phase transition is defined by the behaviour of 
qo(e) = lim; l ^ + q(^,h). For e > we should have qo = 
qbuik, i-e. a small value usually temperature independent; 
at e = we should have qo(0) = Cp and for e < we 
should have <?o( e ) > Gp. This behaviour is present in 
mean field model where metastable states do exist. It 
survives in the real world in the approximation where 
metastable states can be observed. In other words the 
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FIG. 1: A schematic view of the potential W(q) computed 
at various temperatures, decreasing from above to below. 



two replicas system (with the replica a quenched respect 
to the replica t) becomes critical at the same point where 
the dynamics display the mode-coupling singularity. 

One can sharpen the physical picture by introducing 
a potential W(q) defined as follows @, 0, fl8j . We con- 
sider an equilibrium configuration a. We call P a {q) the 
probability that another configuration r has an overlap 
<7o-,t = q- We define 



W(q) — — lim 

v— >» 



v 



(7) 



With probability one when the volume V goes to infinity, 
the potential W(q) does not depend on the reference con- 
figuration a. In other words P a (q) ~ exp(— VW(q)). By 
construction W(qb u ik) = and the vanishing of the po- 
tential W(q) for more than one q value is the distinctive 
characteristic of replica symmetry breaking (this should 
happen below an eventual thermodynamical glass transi- 
tion). In other words we are considering an equilibrium 
configuration a and we define by W(q) the increase in the 
free energy density if we constrain an other equilibrium 
configuration r to stay at overlap q. 

The behaviour of the potential in the mean field ap- 
proximation is described in Fig. [TJ The dynamical tran- 
sition is characterized by the presence of an horizontal 
flex for the potential W (q). Beyond the mean field ap- 
proximation the Maxwell construction should hold and 
the non convex part of the potential disappear, but we 
will not consider this effect. 

A standard assumption is that the behaviour in the 
dynamics mirrors the behaviour in the equilibrium prop- 
erties of two replicas of the system. This assumption is 
usually accepted and there are partial proofs of its va- 
lidity in perturbation theory (we will come later to this 
point by showing how to complete these proofs). Let us 
try to formulate this assumption in a sharp way. 

In the same way as the time may be eliminated para- 
metrically, also the forcing field h may be eliminated in 
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favor of the expectation value of the overlap (i.e. q) and 
all other physical variables can be expressed as function 
of q; this is the usual Legendre transformation of statis- 
tical mechanics. Our assumption is that in the critical 
region (e — > and C —> Cp )the quantities 



We can define two correlations functions 



X4(e|g), T 4 (x -y,e\q) , 



(8) 



computed in perturbation theory in the replica approach 
are the same as those computed in the mode-coupling 
theory (as function of C). This relation holds only in the 
region of q where the forcing field h « 0. 

Let us start the perturbative computation in the replica 
approach. The computations in the mode-coupling ap- 
proach will be rather similar. It is quite clear from Fig. [T] 
that changing the temperature there is a critical point. 
Our task it to compute the critical properties using the 
standard renormalization group tools. We want to de- 
termine the critical exponents 7 and v and the scaling 
function g/Sxj^), that according to the previous assump- 
tion are in the same perturbative universality class in the 
mode-coupling theory and in the replica approach. Let 
us discuss for the moment only the region below the dy- 
namical temperature, i.e. e < 0, since this region is very 
well defined in the perturbative expansion. 

In high dimensions the correlation function is just given 
by the mean field result 0, [H, Hj 



Gi(k,e) 



(| £ | 1 /2 + fc2)2 



thus leading to the following mean field exponents 
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0) 



(10) 



As we will see, there is a crucial difference among the 
double pole behaviour that we find to be valid in the 
general case and the single pole behavior. Notice that the 
double pole form in Eq. (|9t is at variance with the single 
pole behavior predicted in [l9t[2Tj| . if no local conservation 
laws are present. The reasons for this discrepancy will be 
discussed later, when we confirm the presence of a double 
pole by a dynamical analysis. 

In order to compute the critical exponents of the dy- 
namical transition we will start from the perturbation 
theory around the mean field theory. We will consider 
the most divergent diagrams near the critical point and 
we will identify the universality class of the problem. 
In this context we do not see that the critical point is 
in the metastable phase because this would be a non- 
perturbative effect. 

This computation can be done within the replica 
method of 5] assuming a universal cubic effective replica 
action close to the inflection point of W(q) for e — > 0. 
For disordered systems like Potts model, however, this 
model should be in the same universality class of struc- 
tural glasses as far the dynamical transition is concerned. 



Go (a:) = (« o -,T(a:)><g<r > r(0)) - (g<r,r(0)> , (11) 
Gi(x) = (q <7<T (x)qa,T(0)) - (q ( r,T(x))(q (7!T (0)) , 

where the overline denotes the average over the cr's and 
(•) the average over the r's. 

In the case of a fluid with only one type of particles, 
qa.r{x) can be written as <t(x)t{x) (where a(x) and t{x) 
are smeared densities around the point x) and the previ- 
ous equations become 



G (x) = (t(x))(t(Q))o-(x)o-(0) - (t(O)MO) , (12) 
Gi(a;) = (t(x)t(0)) c a(x)o-{0) . 

One finds that near the critical point (neglecting loops) 
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G °( fc > £ ) - ( £ l/2 + fc 2)2 ' Gl ( fc ' £ ) - 



1 



(13) 



In order to compute the loop corrections one must use 
a precise replica setting. In the approach presented here, 
there is an annoying asymmetry among the replica er and 
the replica r. In the replica approach it corresponds to 
take one privileged replica a and n replicas of type r 
and sending the value of n to zero. It is known that a 
partially equivalent formalism consists in using m replicas 
that are constrained to be at overlap q and take the limit 
where m goes to one. In this different formulation we 
have a symmetry Z m in the limit m — > 1 and it may 
be more convenient to use 221. We checked that the 



two formulations are equivalent near the stationary point 
and for simplicity we present the analysis in this second 
formalism. A detailed discussion of this point will be 
presented elsewhere [23J. 

The diagrammatics is the same of an usual <fi 3 theory 
with two differences 12J, |2l|: A) some of the propaga- 
tors have a single pole, others have double poles; B) the 
multiplicity of the diagrams have to be computed in the 
limit where m goes to one and some diagrams give zero 
contribution in this limit. 

One should do a careful analysis: at the end of a long 
analysis 23[ we find that many diagrams give zero contri- 
bution and the dimension where the perturbative correc- 
tions are divergent is 8 [33] . Moreover the diagrams are 
the same of those for lattice animals, with the difference 
that here the effective coup ling constant is positive (for 
lattice animals is negative) [24l426| . The value 8 for the 
upper critical dimension may provide an explanation for 
many of the anomalies found in [27I ] . 

What happens below (and near) 8 dimensions? The 
situation is quite puzzling: the renormalization group 
pushes the coupling constant g 2 toward a large value and 
there is no perturbative fixed point that we can analyze. 
Moreover the terms in the perturbation theory have all 
positive sign and therefore it not easy to estimate the 
sum. This result is not so disturbing. Metastable states 
have finite life time, the free-energy acquires an imaginary 
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part that pushes the singularity in the complex plane; 
in the same way the coupling constant takes an imagi- 
nary part. Although the bare coupling is real, the fixed 
point may correspond to an imaginary coupling constant. 
One could argue that asymptotically the exponent are like 
those of lattice animals for the complex singularity. 

However the previous conclusion may be to hasty. Us- 
ing the same arguments of 24]- 26| one finds that the sum 
of the leading diagrams is related to the solution of the 
stochastic differential equation governing the local fluc- 
tuations of the overlap <f>(x) = q(x) — q (being q the space 
average of q(x)) 

- A^ix) +A + e<p u (x) + g^(x) 2 = u{x) (14) 

where lo(x) is a Gaussian short range noise, that is 
uj(x)ui(y) = 5{x — y). A and g are smooth functions 
of the temperature and they are chosen in such a way to 
implement the condition (j)^ = 0. The two propagators 
are given by the relations 



G (x) = (t>M)M0) > 



d(x) = </> u (xM0) 



-A + e + 2g<Pj x 



, (15) 



where the last equality follows from integration by part 
(and it is correct only in case the solution to the stochastic 
differential equation is unique). 

Neglecting technicalities, the physics is quite clear. The 
choice of the variables a (i.e. the initial conditions in the 
dynamics) induce point dependent shift of the critical 
temperature and the effects of these fluctuations is the 
dominant one. 

One may wonder if there is a direct role of G\ (x) in the 
dynamics. A suggestion is the following. Let us consider 
a theory where the microscopic evolution equations for 
the particle have a stochastic nature. In this case the 
overlap q 0yT (x) will be a function of both the time t and 
of the noise rj(x,t) and it will be denoted by qa- iTI (x,t). 
We can define a different dynamical susceptibility: 



V X 22(t) = (Q CT , n (t) 2 ) - <Q CT , n W) 2 , 



(16) 



and the corresponding correlation G2%{x,t), where the 
overline denotes the average over initial conditions and 
angular brackets average over n 34]. One could argue 
that in the region of time where C(t) is near to (and 
above) the plateaux, G22(x,t) should behaves as G\(x). 
On the other hand in the region where C (t) is small, the 
behaviour of G22(x,t) and G^x.t) should be similar. 

The whole analysis can be redone in the dynamical 
approach where one takes care in an explicit way the de- 
pendence of the correlations on the initial configuration 
(a fact that was neglected in [lj|). The computations 
can be done in a neat way within the Martin-Siggia-Rose 
(MSR) formalism of equilibrium dynamics using a uni- 
versal cubic dynamical action close to the mode-coupling 
transition. Detailed computations shows that taking care 



of the correlations of the initial configuration with the 
evolving configurations one recovers the same result of 
the replica formalism. An explicit isomorphism of the 
two approaches can be shown to be present: the relevant 
computations will be presented in [23] . 

Let us came back to the replica approach and let us 
consider in more detail the stochastic differential equation 
(|LT| that is the resummation of the leading perturbative 
contributions. We arc interested to study it in a non- 
perturbative way. 

It is well know that in perturbation theory the solution 
of the stochastic differential equation is unique and that 
there an hidden supersymmetry [24j-[26J. The supersym- 
metry relates the two propagators and gives 



1 8G (x) 

Gi{x) oc 

X ox 



(17) 



This supersymmetric relation is at the origin of the di- 
mensional reduction: the critical exponents of lattice an- 
imals problem in dimensions D are the same of those of 
the Ising model near the Lee- Yang singularity in dimen- 
sion D — 2. 

However in presence of multiple solutions (as e.g. for 
the Random Field Ising Model) everything becomes more 
complex (multiple solutions cannot be seen in perturba- 
tion theory so that this problem does not affect the per- 
turbative analysis that we have presented above). Super- 
symmetry and dimensional reduction are only valid if we 
average over all the solution with a sign depending on the 
parity of the corresponding Morse index: this weight is 
not the natural one from the physical viewpoint [24l - l26| . 

Which is the correct weight? If we stay within the 
replica formalism it is useful to consider the free energy 



F[ql 



dx 



1 d 2 g(x) 

2 dx 2 



W(q(x))-u(x)q(x)). (18) 



Equation (fl4l can be written as SF[q] UJ /5q(x) = 0. Now 
the natural choice would be to consider among the many 
solutions the one that minimize F[q] u . On the other hand 
this choice is not natural in the dynamics where we would 
like to take the soluti on ci m(x) that maximize q(x). In- 
deed it can be proved [281 ] that there is a solution qAt(x) 
such that qM(x) > q(x) for all x and for any possible 
solutions of the stochastic differential equation (this is a 
well know fact in the Random Field Ising Model (28j|). 
In both cases dimensional reduction is no more valid and 
non-perturbative effects are present. 

A neat formulation of the problem is the following: we 
introduce a fictitious time s and we write the following 
evolution equation: 



dq{x, s) 
ds 



-Aq(x,s) + W'(q(x,s))-uj{x) (19) 



with the boundary conditions at time s = given by 
q(x,0) = 1. In the region below T c , where the equation 
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(TTJ| has many solutions, the solution relevant for the dy- 
namics is uniquely identified as 



q*{x) = lim q u (x,a) 



(20) 



We are near to the end of our journey. We still have to 
compute the critical exponent and the critical behaviour 
of the correlation of (x) . Techniques introduced in [2^] 
could be used to achieve this goal. 

We remark that the fictitious time s in the equation 
(fHl]| should not be identified with the real time: the de- 
pendence of q(s) if we solve equation (1191) is quite different 
from the behaviour of q(t) in mode-coupling theory and 
there should be no confusion among the two variables. 
However one may make en passant the conjecture that, if 
we consider only reparametrization invariant quantities, 
the mode-coupling equations and eq. (fT9l) stay in the 
same universality class: in other words, if in the mode- 
coupling theory we eliminate the real time t in favour 
of q and we write everything as function of q and in eq. 
(I19[) we eliminate the fictitious time s in favour of q, the 
dependence of the physical quantities on q should be the 
same near the critical point for the models. This con- 
jecture is trivially valid if there is only one solution of 
equation (fT4")) . The point is to understand its correctness 
beyond perturbation theory 

This conjecture may be generalized by introducing a 
more general evolution equation that can be used to com- 
pute the properties of reparametrization invariant quan- 
tities when activated processes are present below T c : 



dq(x,s) 
ds 



-Aq(x, s) + W'(q(x, s))-cj(x)+T](x, s) (21) 



where n(x, s) is a thermal noise. The discussions of the 
consequences of these conjectures cannot be done here, 
however it has not escaped to how attention that eq . (Effl) 
can be used to explain the experimental results of [30j|. 

Summarizing, we have clarified the predictions and the 
limitations of the perturbative expansions for the critical 
properties of glasses, finding an explicit mapping among 
the replica formalism and the mode-coupling approach 
in the framework of the MSR approach to the dynamics. 
We have also conjectured a mapping with a differential 
stochastic equations that should be valid beyond pertur- 
bation theory, whose interesting consequences should be 
studied in details. 
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